Exact Solution of a Yang-Baxter Spin-1/2 Chain Model and Quantum Entanglement 
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Entanglement is believed to be crucial in macroscopic physical systems for understanding the 
collective quantum phenomena such as quantum phase transitions. We start from and solve exactly 
a novel Yang-Baxter spin-1/2 chain model with inhomogeneous and anisotropic short-range inter- 
actions. For the ground state, we show the behavior of neighboring entanglement in the parameter 
space and find that the inhomogeneous coupling strengths afi^ect entanglement in a distinctive way 
from the homogeneous case, but this would not affect the coincidence between entanglement and 
quantum criticality. 
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I. INTRODUCTION 

Entanglement was conventionally considered to be a 
quirk of microscopic objects and has been recognized to 
be ubiquitous and so robust that it promises applica- 
tions in the quantum communication and computation 
like technologies [11]. In the last few years, there was an 
increasing interest in entanglement in macroscopic phys- 
ical systems 0, Q- Entanglement may lead to further 
insight into condensed matter physics. For example, in 
statistical mechanics, given that quantum phase transi- 
tions (QPTs) occur at absolute zero and are driven by 
quantum fluctuations, entanglement may provide addi- 
tional correlations for QPTs [J, [^] that have no classical 
counterpart. In return, materials and experience built up 
over the years in condensed matter are helping in finding 
new protocols for quantum computation and communi- 
cation. 

Studying entanglement of the ground state in macro- 
scopic physical systems is crucial to understand a large 
variety of collective quantum phenomena. For a one- 
dimensional spin-1/2 XY chain Q with short-range in- 
teractions 
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(1) 
the entangled degree (ED) between any two nearest- 
neighbor particles keeps the same for the translational 
symmetry, and its derivative is capable to fulfill the role 
of an order parameter to characterize QPT at the critical 
point A = 1 d, [^, 0] ■ Potential as it is, we should ask 
whether such an observation is universal enough to assure 
all correspondences between entanglement and QPT. For 
example, in Ref. [8] it demonstrates a long-distance en- 
tanglement appearing for values of the microscopic pa- 
rameters which do not coincide with known quantum 
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critical points. In addition, if the short-range interac- 
tions are not homogeneous, e.g., a dual chain, ED be- 
tween two nearest-neighbor particles would generally not 
keep in accordance for different sets of chains. At this 
case, questions arise on what ED is and whether it well 
coincides with critical points. 

For self-contained and later use, let us recall the Yang- 
Baxter approach for entangled states with two qubits 
9]. The unitary R(9,(j)) matrix in this approach, tak- 
ing R{9) ~ sin + cos OM {(f>) with the generalized four- 
dimensional imaginary unit A4 {Ai^ = —1), is used to 
produce entangled states when acting on direct-product 
states and EDs of the resulting states are simply | sin 20 1. 
Supposing that 6 is time-independent and (j) is time- 
dependent, we can get a Hamiltonian through H{0, 0) = 
ih{dR/dt)W (see Eq. ([2])), which only governs the evo- 
lution of entanglement varying with parameter 9. At the 
case oi 9 = 7r/2, R corresponds to an identity operation 
with no entanglement at all and H thus vanishes. Such 
a peculiar behavior of H actually creates a nonanalytical 
point of the ground state energy with respect to 9 and 
this should be reflected by some properties of its ground 
state such as the geometric phase (GP). When extended 
to an infinite lattice, the possibilities are richer and the 
vanishing point may correspond to a QPT. 

The purpose of this paper is twofold: one is that 
we solve exactly a novel Yang-Baxter spin-1/2 chain 
model with alternating coupling strengths by means of 
the Jordan- Wigner transformation and GP of the ground 
state is examined for QPT; the other concerns the con- 
sequence of entanglement between two local nearest- 
neighbor particles of the chain, based on which we check 
whether entanglement under inhomogeneous coupling 
strengths can well characterize critical phenomena. This 
article is organized as follows. In Sec. |lTl we introduce 
and exactly solve an inhomogeneous Yang-Baxter spin- 
1/2 chain model. Based on the solution, we investigate 
the quantum criticality by analyzing GP of the ground 
state, and study the effect of inhomogeneity on entangle- 
ment between different nearest-neighbor sites and next- 
nearest-neighbor sites in Sec. IIIIl At last. Sec. IIVI is 



dedicated to the conclusion. 



II. THE YANG-BAXTER SPIN-1/2 CHAIN 
MODEL AND EXACT DIAGONALIZATION 

The Yang-Baxter equation was originated from solv- 
ing the (5-function interaction model by Yang [10| and 
the statistical models by Baxter [ll|, and was then in- 
troduced to solve many quantum integrable models by 
Faddeev and Leningrad Scholars ^^. It plays a funda- 
mental role in the theories of 1-1-1 and 2-1-1 dimensional 
integrable quantum systems, including lattice statistical 
models and nonlinear field theory. For example, Yang's 
^-matrix in YBE for the nth and {n + l)th particles is 
R{u)nji+i = l + u'Pn,n+i (u is a spcctral parameter, i.e., 
one-dimensional momentum and V is permutation satis- 
fying V^ = 1); it yields the XXX chain model through 
Hn,n+i ~ dR/ du\u=Q, when V takes its four-dimensional 




representation of Vn,n+i 



[1 + cr n- (7 n+i) with a being 



the Pauli matrix. Considering the form of R{On) in Sec. 
m it has a two-body interacting Hamiltonian [9|] : 

Hn,n+1 = -fv^ cos On [ COS 6'„ [S^ -f S'fj+i) 

+ su,en{e^StS+^, + e-^^S-S-^,)\, (2) 

where S^ ^ S"" ± iS^ with S^'?''^ = cr^'2''V2, and cj) 
is the flux dependent of time t and it takes (j){t) = tut, 
denoting procession angle of spins around the z direc- 
tion in a rotating magnetic field. For many particles, we 
should sum all of them as H — ^^^ Hn,n+i and if all of 
9n are taken to be the same, it would correspond to a 
homogeneous chain, otherwise it would correspond to an 
inhomogeneous one. 

From Eq. ([2]), the family of Hamiltonians that is pa- 
rameterized by (f) is clearly isospectral, and, therefore, the 
critical behavior is independent of cj). In fact, we can see 
that the spin raising-raising or lowering-lowering struc- 
ture in Eq. ^ allows a rotation for each spin around 
z-axis and such a rotation transformation can be em- 
ployed to adjust the value of phase factors in Eq. ^, e.g., 
n - g{^/2)Hg^i^/2) and 5(0) = UZi^'"'''^^' giving 
g{H2)S+g\<l>/2) = e~^^/^S+ and g{<l>/2)S-g^W2) = 
e^'^'^S„ . Then the Hamiltonian is reduced to 



n = --fiw^cos6l„[sin6'„(«+i-cr>f^+i) 

n 

-Kcos6'„«-}-<+i)]. 

Comparing H with the XY chain in Eq. ([T|), one sees 
that the structure of dominant two-body interaction in 
Ti. is exactly that in Hxy yet under 7, A ^ 1 limit. Like 
HxY, there is a global Z2 symmetry for H that keeps 
it invariant under a unitary transformation Hn'-'^n- In 
the following, we will see that the ground state does not 
break such a symmetry. 

By noting the fact that it can define these opera- 
tors: J^ = 2\^n + '^n+l)' Jn — 2\^n ^n+l + ^n ^n+l) 



FIG. 1; A chain model with two-body interactions having the 
Hamiltonian equation Q : the parameters 6\ on all solid dark 
lines and 62 on all dashed red lines describe two different cou- 
pling strengths; each ellipse including solid and dashed ones 
represents a two particle composite qubit with the angular 
momentum Ji. 



and jy = jiS^S^^i — S'^5'^+i) that satisfy the an- 
gular momentum commutation relation, we see that 
Jn actually only occupies a subspace spanned by | tt 
)„,„+! and I ii)„,„+i that belongs to a j = 1/2 angu- 
lar momentum representation. Thus we can write the 
Hamiltonian into an effective NMR-like form Hn,n+i = 
—Bn{t) ■ Jn, where the magnetic field is Bn{t) = 
2hujcos6n{sm6n cos(/)(t), — sin6'„ sin (j>{t) , cos 6 n) ■ Its 
eigenvalues are readily given by £"„ „^]^ = ±Huj cos 9n 
and eigenstates in accordance are \E^) = cos^| tt 
) ± sin^e-'-^l ii), both EDs of which are |sin6l„|. In 
the vicinity of 6'„ — 7r/2, we can set cos6'„ = S and see 
that it keeps the forms of the eigenstates \E^{S)) and 
hence EDs even when ^ ^ 0. 

Now, let us extend Eq. ^ to an inhomogeneous chain 
which approaches infinite when the particle number takes 
to be arbitrarily large. For simplicity, we confine bonds 
between any pair of odd-even numbered nearest-neighbor 
sites to be the same and characterized by parameter 9i 
and that between any pair of even-odd numbered nearest- 
neighbor sites to be the same and characterized by 6*2 (see 
Fig. [T]). By requiring the periodical boundary condition, 
there are totally 2N sites and the Hamiltonian can be 
written as 



N 

H = -Y,{Bi-J2n- 

n=l 



-1 + ^2 • Jin), 



(3) 



which can be interpreted simply as two sets of spins in 
an external field with different coupling strengths along 
z-axis. The composite qubit Ji satisfies [J^, Jj] — for 
\i—j\ > 2. At the same time, the property of H — zero in- 
teraction at zero field — provides a critical phenomenon at 
Oi = 02 = Ti 12 and is similar but not exactly the same to 
the on-site exchange interactions ^ISj and superexchange 
interactions [ij with ultracold atoms in optical lattices, 
which would vanish as the optical field is taken off. Here 
we will focus more on the entangled aspect of these in- 
teractions by using the Hamiltonian Eq. ([3]), which has 
a closed relation with entangled states. As for this, it 
is interesting to ask what ED is when extending to an 
infinite chain, especially for an inhomogeneous chain and 
whether it can effectively determine the critical point. 

To solve the model, let us first introduce the follow- 
ing Jordan- Wigner transformation to represent spin op- 
erators at sites with spinless fermion operators: a„ = 



(]!,<„ erf )S'+ or S+ = a„e-"^'<"°'°'. Then, note that 
the summation in Eq. ([3]) is either on aU even or on aU 
odd indices, and so we should distinguish even and odd 
to define their corresponding forms in momentum space, 
respectively, as 



same symmetry as the Hamiltonian. Alternatively, the 
ground state can also be expressed by 



\9) - n {^ke'^^'\0)-km + «.e-^/2|l)l 



fc|l>fe 



(8) 



1 ^ 

m— 1 



1 



N 



N ^ 

m— 1 



«2m-l, 



(4) 
where the reduced momentum k = —M, . . . ,M with 
M = {N — l)/2 for TV odd and fermion operators 






(t) „o(t) 



) anticommute with each other. Thus the Eq. 



([3]) can be written into 

H = -^hw J2 [iCke^'^alaU + h.c.) 

-A{al^al + al^al-l)], (5) 



M 



with ^k = sin 2826'^'''''/^ - sin 26*1 and A = cos 26*1 + 
cos 202 + 2. The Hamiltonian H can be diagonalized by 
using the Bogoliubov transformation and the result is 



k 



(6) 



The eigenspectra contai n two ba nds of quasiparticle 
excitations: e^ = ±-\/|$fcp + A^. The transformed 
fermion operators a^ — UkC^'^^^a'j^ + Wfee^"^/^a'l\. and 
(3k = Uke'">'/'^al + Vke~''''>'/'^a'\, where Uk = -Uk = 
(A + e±)/[2e±(A + e±)]i/2 and Vk = vl = -a/[24(A + 
e^)]-'^/^ for different bands e^. For these coefficients, 
there is m^m^ + Vkvl = 0. 

The ground state \g) of H is the vacuum of the 
fermionic modes, satisfying ak\g) — and (3k\g) — for 
all fc. Generally it is hard to write the ground state ob- 
viously into a spin superposed state, but for our model 
there is a fact that the Ned state— |^i) = | ]l)®^ 
or 1^^) = I It)®^ — just corresponds to the zero en- 
ergy eigenstate of H. Take |^^) for example: it has 
ak\^^) = and therefore it is identical to |$^) = 
nfe/^fels)' which inversely gives an expression for the 
ground state \g) = Hfe /3fc|*^), i-e.. 



b)-n{E 



e N ' 



k 771 — 1 
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Vke 
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( n '^n^2™-i 



(n <^f)sL]}\n) 



l<2m 



l<2m-l 



(7) 



which would return to the biparticle case (i.e., \E )) 
if one takes iV = 1 and 6*2 = 7r/2. When 6*1 ^ 7r/2 
and 62 = 7r/2, we can see that Eq. ^ becomes to 
describe N isolated dimers, which has an exact ground 

state, nrr7(cOS^| TT)2r7i-l,277i+sin^e-"^| ii)2m-l,27n)- 

From Eq. ([7]) , it can be seen the ground state is invariant 
under the global Z2 transformation and so it keeps the 



where lO)^'*^ and \^Yk'^ are the vacuum and single exci- 
tation of the fcth mode, a^'*^, respectively. The ground 
state is a tensor product of states, each lying in the 
two-dimensional Hilbert space spanned by |0)° j,|0)| and 
|l)-fe|l)fc- Such a form provides us a convenient way to 
discuss its dynamical property such as GP of the ground 
state. 



III. QUANTUM CRITICALITY AND 

ENTANGLEMENT 

QPT occurs at a point in the external parameter space, 
where there can be a level-crossing and excited levels be- 
come the ground state, creating a point of nonanalyticity 
of the ground state energy as a function of external pa- 
rameters jT^. With the Hamiltonian Eq. ([6]) in consid- 
eration, we take 9i^2 as those external parameters. Ob- 
viously, at the point of 6*1 = 6*2 = 7r/2 all energy levels 
cross and hence it is a critical point, but it is different 
from the conventional QPT by having a vanishing Hamil- 
tonian and for convenience of discussion, we might as well 
caU it a QPT. Recently, GP of the ground state [ll[l3 
and ED between nearest-neighbor particles [J,la| for a ho- 
mogeneous Heisenberg XY chain were proposed to char- 
acterize the criticality of QPT. In this section, we inves- 
tigate GP and ED for our novel inhomogeneous chain, 
analyze their behaviors as the parameters 61^2 vary, and 
further discuss their nonanalytical property in the prox- 
imity of QPT. 



Quantum Criticality Characterized by 
Geometric Phase 



GP of the ground state, accumulated by varying 
the angle (p from to 27r, is described by /3g — 
W lo^ i9\'^'^'t>\9)'^'f'j and by utilizing Eq. ([5]) it is 



fl 



^Bi 



N 



Uk\ 



m 



^E^/^t 



N 



(9) 



with /?+ - 



g P~ . The term Pk = """^/^fc is a geomet- 

ric phase for the fcth mode, and represents the area in 
the parameter space enclosed by the loop determined by 
(6*1,02,^)- One can see that when we turn off the cou- 
pling between dimers by setting 62 = 7r/2, GP would re- 
turn to the biparticle dimer case with f3g — 7r(l — cos 9i). 
To study quantum criticality, we are interested in the 
thermodynamic limit when the spin lattice number N — > 
00. In this case the summation jj J2k=-M '^^^ ^^ ^^~ 
placed by the integral i J^ dip with ip — 2nk/N; GP in 




concurrence for the state of the ith and jth particles is 
defined as 



FIG. 2: (Color online) The diagram for the Berry phase of 
the ground state of the spin-1/2 chain model: (a) the Berry 
phase /3^ corresponding to the ground state; (b) its deriva- 



tive dPg /d92 as a function of 9\ and 
8^(3^/861862 as a function of Oi and 62- 



the thermodynamic limit is given by 



its derivative 



Jo 



(10) 



where the energy spectra e^ ~ ±-y|^^^p + A^ with 
|^^|2 ^ sin^ 26*1 + sin^ 26*2 - 2sin26'i sin 26*2 cos (p. 

To see the quantum criticality obviously, we plot GP 
l3g and its derivatives 8(5/ 862^ 8"^ (i / 8B\dQi in the param- 
eter (Q\,B-2) space, shown in Fig. [2l It can be seen, from 
Fig. IHa), that there is a conical intersection at the point 
Q\ = Qi = T^ I'i'i which indicates a nonanalytical point 
there. The nonanalytical property at the critical point 
can be seen obviously from the diagram of GP derivative 
[see Fig. [2Ic)]. However, as we pointed above, such a 
QPT point is trivial, since at the point the whole Hamil- 
tonian vanishes and hence it appears to be exotic there. 
If fixing one parameter, say 0i, it would correspond to the 
uniparameter case and we should check whether there are 
other critical phenomena as varying Qi- The derivative 
of GP with 02 is plotted in Fig. [U^b), from which we can 
see, except for the above critical point and its vicinity, it 
is analytic everywhere. So there is no additional critical 
point. 



B. Entanglement of the Ground State 

In this section, we confine our interest at entanglement 
between nearest-neighbor sites in the chain, given long- 
distance entanglement decays rapidly with the distance 
(see below). To describe entanglement, we use the con- 
currence [Ijl of a biparticle state, related to the "entan- 
glement of formation" [19], to define ED of a state. The 



C(i,j) = max{ri -r^-r^- r4,0}. 



(11) 



where ri_2,3,4 are the square roots of the eigenvalues of 
the product matrix R = p{i,j)p{i,j) in descending or- 
der; p{i,j) is the density matrix of the ith and jth spin- 
1/2 particles and the spin flipped matrix is defined as 
p{i,j) = (T^ (8) (7^ p*{i,j)(J^ (81 cr^. If it is a pure state, e.g. 
\E^), the density matrix p{n, n+l) = |_E^)(£'~| and the 
concurrence quantifying entanglement is C{n,n -I- 1) = 
I sin 6n\. If it is a biparticle state in a multiparticle sys- 
tem, p{i,j) would represent a biparticle mixed state re- 
duced from the multiparticle density matrix p. 

For the chain in consideration, translation invariance 
of dual lattices implies that C(2m, 2m + 1) = Ce(l), 
C(2m - 1, 2m) = Co(l) and C(2m - 1, 2m + 1) = C(2) 
for all m. The concurrence will be evaluated as a func- 
tion of the relative position |« — j| between the ith and 
jth spins and parameters 0i^2- All information needed is 
contained in the reduced density matrix p{i,j) obtained 
from the ground-state wavefunction after all the spins 
except those at positions i and j have been traced out. 
The resulting p{i,j) represents a mixed state of a biparti- 
cle system. The structure of the reduced density matrix 
is obtained by exploiting symmetries of the chain. The 
nonzero entries of p{i,j) can then be related to the vari- 
ous correlation functions [20, 0, u3 as 



pihj) = (^+K)'^f '»i + K)i®^i + ('^r<^i) 



X,Y=x,y 






W®a])/A. 



(12) 



Correlation functions under the ground state can be 
evaluated by using the fermionic representation and the 
simple identity 1 — 2a,Jja„ = (a\ -\- an){a\ — an)- For each 
pair of fermion operators {a^a^)^ we can further define 
two Majorana fermion operators (A, J5): A„ = a\ + an 
and Bn — i{a\ — an) with A^ — A and B'^ — B. 
Exploring them, we can write the Pauli matrices as: 
< = Y{i<n[i-^)AiBi]An, al = ni<n[(70^iSi]Sn, and 
CT,^ — {~-i)AnBn- A two-body correlation function, say 
('''m'''n)|m<nj Under the ground state, is 

(«)=~*(s™ il ^-^MBi)a\ (13) 

= ((-i_B„JA„+i(-i_B,„+i) • • ■ An-\{~iBn-\)An). 

Since the expectation values are with respect to a free 
Fermi theory, the expression on the right-hand side can 
be evaluated by the Wick's theorem 15, 23], which re- 
lates it to a sum over products of expectation values of 
pairs of operators, i.e., (AiA„J, (B;B,„), and {BiA.m). 
The evaluation of average values of these pairs is dis- 
played in Appendix \^ In order to see the macroscopic 
property of entanglement, we define two fc-independent 



functions: 



k 

G{m-n) = l^e-^'*^("-™)2wfcz;fc, (14) 



which have summed all frequencies in the momentum 
space to be the form in the position representation. We 
can see, J^{(i) is nothing but the one proportional to 
GP in Eq. ([9]). In this respect, we may well say T^Q 
are macroscopical quantities, which under the thermody- 
namical limit N —f oo can be calculated still by making 
the replacement jjj^k^'^lo ^'^■ 

Next, we would focus on entanglements between the 
nearest-neighbor sites. In principle, the numerical results 
on the concurrence of such two-site density state can be 
performed readily according to its definition introduced 
above, but before that, we find the concurrence depends 
only on the above two fc- independent functions J^{\n — 
TO I) and Q(rn — n). As an illustration, we give out the 
expressions of ED in the form of concurrence between 
odd-even and even-odd neighboring sites, respectively, by 

Coil) = max{0, 1^(0)1 -i|.F(0)2 + g(0)2-l|}, 

Ce(l) = max{0,|t;(l)|-l|.F(Of + g(lf-l|}.(15) 

Their varying trends in the parameter space are displayed 
in Fig. [nja) associated with contours in EJb). Note 
that although we plot merely the concurrence of even- 
odd neighboring sites in Fig. [3l that of odd-even neigh- 
boring sites can be immediately gotten by noticing the 
symmetry: Co{Oi <-^ 6*2) = Cg. 

Fig. Elja) shows that the maximum of Ce(l) (or Co(l)) 
approaches but cannot reach one (i.e, the intrinsic maxi- 
mal value of ED) because of the vanishing Hamiltonian at 
the value. At 9i = 7r/2 even-odd pairs become decoupled 
from the rest and the behavior of Ce for 62 — » 7r/2 is the 
same as in the biparticle case. Whereas, the difference 
refers to the minimum of Ce(l) that relates to disentan- 
gled states with vanishing ED and from Fig. Efb) we can 
see that in the parameter space it occupies a considerable 
regime enclosed by the dashed red lines for disentangled 
states. Thereby, at the interface of the entangled and 
disentangled states it will correspond to the unsmoothcd 
part of ED mainly due to the unsmoothed definition of 
the concurrence [see Eg. pT]) ] and the inhomogeneous in- 
teraction strength in the chain, which would as well im- 
ply the nonanalytical behavior for derivatives of ED with 
respect to 6*1.2. In contrast, at the homogeneous case of 
Ce — Co achieved along the diagonal line of cither 9i — 62 
or 6*1 -1-02 = TT in Fig. El^b), ED has only one unsoomthed 
point which corresponds to the critical point appropri- 
ately. 

Fig. ^c) and (d) show respectively the first and second 
derivatives of ED with respect to parameters 9i^2- The 
dominating divergence at the central region is well in 




FIG. 3: (Color online) The diagram for concurrence in the 
parameter space: (a) ED in the form of concurrence is shown 
and Ce(l) denotes entanglement between even-odd numbered 
two nearest-neighbor sites; Its corresponding contours are 
also shown in (b); (c) and (d) display the first derivative 
dCe{l)/d92 and the second derivative d^Ce{l)/ 861862 of 
Ce(l), respectively. 



agreement to that appeared in the derivatives of GP and 
thus indicates QPT. Except this, it does not diverge at 
other points, exhibiting instead just a finite discontinuity 
at the border between separable and pairwise entangled 
sectors. This type of non-analytic behavior is quite dis- 
tinct from that at the critical point and stems just from 
the definition of concurrence. In addition, the divergence 
includes both upward and downward directions and this 
reflects the increasing and decreasing trends of ED in the 
vicinity of the nonanalytical region. 

As for two next-nearest-neighbor sites and other far- 
ther neighboring sites, ED should decay rapidly with 
the distance (generally even more rapidly than standard 
correlations) for the short-ranged interaction. To illus- 
trate it, our numerical result demonstrates that the next- 
nearest-neighbor concurrence C(2) vanishes for the whole 
chain. As a result, it is sufficient for us to only take the 
nearest-neighbor biparticle entanglement into account on 
this chain, while for multiparticle entanglement it may 
has a connection with the so called topological quantum 
phase transition (e.g., [2J|) but not discussed here. 



IV. CONCLUSION AND DISCUSSION 

To summarize, we have demonstrated an exact solu- 
tion to a particular spin- 1/2 chain model with alternat- 
ing nearest-neighbor coupling strengths and have ana- 
lyzed the influence of inhomogeneous interaction on the 
ground state through GP and ED approaches. By eval- 
uating GP of the ground state, we display its behavior 
at the parameter space, from which a critical point could 



be determined through the divergent derivative of GP. 
After that, via examining the biparticle entanglement 
by virtue of concurrence, we also show the tendency of 
ED with respect to parameters and find that ED and 
its derivatives can determine the critical point as well as 
GP does at the inhomogeneous case. Although ED has 
an unsmoothed definition at the border of separable and 
pairwise entangled sectors, it does not exactly affect the 
nearest- neighbor concurrence to detect the critical point, 
exhibiting instead a clear signature of it in its derivative 
as seen in Fig. [2{c). 

As remarked earlier, the specific parametrization of the 
Hamiltonian has an intimate relation with ED at the bi- 
particle case and when extending it to an infinite lat- 
tice, the result is interesting: for homogeneous coupling 
strengths, ED has suppressed values with its maximum 
far less than one, keeping equal for every pair of nearest- 
neighbor sites; for inhomogeneous coupling strengths, ED 
appears to have different values between even-odd num- 
bered and odd-even numbered nearest-neighbor sites, 
from which a very high ED is available (see Fig. ^. 
In a way, the property might apply to other inhomoge- 
neous lattice models as a manifestation of general princi- 
ples. Also, the analysis of the inhomogeneous entangle- 
ment for a condensed matter system is possibly of great 
importance for creating ideal entanglement resources in 
quantum information processing. 
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APPENDIX A: CALCULATION OF 

EXPECTATION VALUES OF PAIRS OF 

OPERATORS 

In determining the expectation values of pairs of Majo- 
rana operators under the ground state, we have used the 
representation (afc,/3fc) and the definition of the ground 
state ak\g) = Pk\g) = 0. If we define two independent 
real functions, i.e., Eq. ([14]) by 

:F{\n~m\) = l^e^^("-™)(|u,.p-|f;fcn 



fe 

= ]YXl'f^^'^2^2COs[— /^(TO-r^-l)]- 
fc 
sin20l cos[— fc(m — n)]}/€i^ , 

which are obtained under the thermodynamical limit 
A^ ^ oo by making the replacement -^ X)fe ~^ ^ /J^ '^'^' 
then the average values of Majorana operator pairs can 
be written as 

{B2nB2m) = (-B2ri-li?2m-l) = Sn,m, 

{A2n{-iB2„i)) = {A2n-l{-iB2,n-l)) = ^{\n - m|), 

{A2nA2m-l) = {{-iB2n){~iB2rn-l)) = i siu (/)t/ (m - n) , 

{A2n{-iB2rn-l)) = {{-iB2n)A2m-l) = COS (jiQ (m - u) . 

From these average values, we can calculate all correla- 
tion functions by using Wick theorem. For the nearest- 
neighbor case of p{2m — 1, 2to), we have 

Km-l<^fm> = ((-*^2m-l)A2™) = - cos 05(0), 



'2m-l'^2m/ 



.V \ - 



'{A2rn-l{-iB2ra)) = COS0t/(O), 



'2m-l'^2n 



.V \ - 



'2m-l"2ml 



sin(/.6;(0), 



J^ Yl cos[^/c(n - m)] A/e±, 



{^2m-1^2m) — (^2m-l ( — *-S2„i-l)^2m (" *^2m)) 

For the nearest-neighbor case of p{2m, 2m + 1), we have 

{'^2m'^2m+l) = ((- JS2™) Aa^+l) = COS (/)5(1), 

WmC^L+l) = -{A2,ni-iB2„,+l)) ^ ~COS(l>g(l), 

Kmf^L+l) = «„CT2m+l> =- sin 05(1), 

(^2m^2m+l) = (^2m (-i-B2,„) ^2™+! (-ii32m+l)) 

= .F(0)2+ 5(1)2. 

At last, for the next-nearest-neighbor case of p{2in — 
1, 2m -1-2), we have 

{'^2m-l'^2m+l) = ((-*52m-l)^2m(-«S2m)^2m+l) 

= -5(0)5(1) -.F(0).F(1), 

i'^2m-l'^2m+l) = ^ (^2m-1^2m (^ *-S2m) (" *-S2m+l)) 

= -5(0)5(1) -.F(0).F(1), 

\'^2m-l'^2m+l/ = \'''2m-l'''2m+l/ = '-'' 

W2m-l'^2m+l) = (^2m-l ( — «-B2,„-l) A2m+1 ( — i-B2m+l )) 

= .F(G)2 - .F(l)2. 
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